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Dark matter charged under a new gauge sector, as motivated by recent data, suggests a rich 
GeV-scale "dark sector" weakly coupled to the Standard Model by gauge kinetic mixing. The new 
gauge bosons can decay to Standard Model leptons, but this mode is suppressed if decays into lighter 
"dark sector" particles are kinematically allowed. These particles in turn typically have macroscopic 
decay lifetimes that are constrained by two classes of experiments, which we discuss. Lifetimes of 
10 cm < CT < 10^ cm are constrained by existing terrestrial beam-dump experiments. If, in 
addition, dark matter captured in the Sun (or Earth) annihilates into these particles, lifetimes up to 
~ 10^^ cm are constrained by solar observations. These bounds span fourteen orders of magnitude 
in lifetime, but they are not exhaustive. Accordingly, we identify promising new directions for 
experiments including searches for displaced di-muons in B-factories, studies at high-energy and 
-intensity proton beam dumps, precision gamma-ray and electronic measurements of the Sun, and 
milli-charge searches re- analyzed in this new context. 



I. LIGHT LONG-LIVED PARTICLES FROM 
NEW GAUGE SECTORS 

Recent astrophysical anomalies have motivated models 
[H EJ [3l m of TeV-scale dark matter interacting with 
new states at the GeV scale. A compelling candidate is a 
new gauge force, mediated by a GeV-mass vector-boson 
that mixes kinetically with hyper charge: 



C 



1 



(1) 



4 



2 cos 



w 



-'Dark- 



A program of searches in high-luminosity e^e~ collider 
data m [71 18] , new fixed-target experiments [9l [TOj , and 
high-energy colliders [TTJ [12] can search quite exhaus- 
tively for the new gauge boson A' if it decays into a lepton 
pair through the e-suppressed kinetic mixing term. 

The A' can, however, be accompanied by "dark sector" 
particles (other vectors, Higgs-like scalars, pseudoscalars, 
or fermions) to which it couples without e-suppression, 
so that the decays of A' into dark sector particles domi- 
nate. These particles may decay to Standard Model mat- 
ter (and indeed some must^ if dark matter annihilation 
is to explain the astrophysical anomalies) , but unlike the 
A' their lifetimes are typically macroscopic — we will call 
such particles "long-lived particles" (LLP's). 

In this paper, we discuss observational constraints on 
the decays of the LLPs which cover lifetimes spanning 
fourteen orders of magnitude, 10 cm ^ cr < 10^^ cm. A 
schematic diagram illustrating the type of experiments 
we consider is shown in Fig. [T] The constraints fall into 
two categories: 

Terrestrial beam dump experiments (Sec. [ll|: Ex- 
periments at proton and electron beam dumps 
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FIG. 1: A pictorial summary of the two similar approaches 
discussed in this paper: A strong source, either a beam dump 
or dark matter annihilation in the Sun, produces dark sector 
particles. LLP's can penetrate the shield, either the Sun's 
interior or the beam dump shielding. Decays downstream 
over a baseline, ^Baseline' detected. 



search for decays of long-lived particles in an in- 
strumented region downstream of a thick shield. 
Though A' production scales as e^, the integrated 
luminosities ~ 10^ ab~^ allow searches for LLP's 
with mass ^ few GeV and 10 < cr < 10^ cm. 

Annihilation of Dark Matter in the Sun (Sec. 



Ill): Dark matter captured in the Sun through 
A^mediated scattering annihilates efficiently, with 
a rate scaling as e^. Annihilations into LLP's 
lighter than the r threshold result in e, /i, and 
light hadrons, which can be detected when the 
LLP decays outside the Sun but before passing the 
Earth. In this case, neutrino, electron, and gamma 
ray observations together constrain dark matter 
annihilation into LLP's with 10^ cm ^ cr < 10^^ 
cm. 



In Sections [TT| a nd |III[ we describe the best available lim- 
its. In Section |IV[ we interpret these limits in a par- 
ticular illustrative case: a Higgs-like scalar in the dark 
sector that decays into two Standard Model fermions. 
We conclude in Section |V| and highlight the value of 
direct A' £~^£~ searches, searches for displaced di- 
muons in B-factories, high-energy and high-intensity pro- 
ton beam dumps, precision gamma-ray and electronic 
measurements of the Sun, and milli-charge searches in 
the context of these models. 
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FIG. 2: Feynman diagrams for (a) A' decay into a dark-sector 
higgs and (b) higgs^-strahlung process. 



A. Decay Lengths in Dark Sectors 

Before proceeding, we summarize the parametric scal- 
ing that governs the decays of various particles in a dark 
sector which is coupled to the Standard Model through 
gauge kinetic mixing. 

In non-Abelian dark sectors, multiple gauge bosons can 
mix with one another [H [121, so that each acquires an 
effective kinetic mixing eg// and lifetime oc ^^ff 
tides X of other spins can have three-body decays into 
Standard Model charged matter (x f~x' through 

an off-shell A')^ with lifetimes cx but longer than 
those of vectors due to three-body phase space. For 
6 ~ 10~^ — few X 10~^ and moderate mass hierarchies 
these lifetimes are less than 1 m, and multi-body final 
states can be observed in a collider detector. B-factory 
searches in higher- multiplicity final states (e.g. [13]) can 
be quite sensitive to these decay modes [6^, 7^ . 

Much longer lifetimes can arise from scaling or de- 
cays through higher-dimension operators — these are the 
modes of greatest interest for our discussion. In partic- 
ular, we consider a Higgs-like boson ho below the A' 
mass. The A' can decay into Hd only in models with 
multiple dark- sector Higgses, as the second decay prod- 
uct must be a pseudoscalar. Several different paramet- 
ric decay lifetimes are possible. The dark Higgs can de- 
cay through gauge kinetic mixing, with scaling, or 
faster if it mixes with Standard Model Higgs bosons (e.g. 
through a quartic term ex\hD\^\h?\)- If the Dark Higgs 
is stabilized by an accidental symmetry so that its decay 
is mediated by higher-dimension operators as noted in 
[14l[^, the lifetimes can be significantly longer. 

For suitable parameter choices, all three decay scenar- 
ios (kinetic mixing, higgs mixing, and higher-dimension 
operators) lead to lifetimes 10 cm ^ cr < 10^^ cm, 
so that beam-dump and solar limits apply. The ki- 
netic mixing decay is a simple and illustrative bench- 
mark, because the decay lifetime and terrestrial produc- 
tion cross-sections both depend on e and the particle 
masses, with no additional free parameters. After dis- 
cussing the model-independent limits, we will consider 
this case explicitly in Section [iVj 



II. BEAM-DUMP EXPERIMENTS 

Production of dark-sector states in proton beam-dump 
experiments was discussed in [10], which identified three 
production modes (meson decays, "higgs'-strahlung" , 
and resonant A' production). For concreteness, we will 
assume that a dark sector scalar higgs /i/) is the LLP 
of interest but the results in this section are insensitive 
to details of the model, except for e characterizing the 
production cross-section, the A' mass and the mass and 
lifetime cr of the LLP. We will consider two reactions 
for producing scalar dark sector particles: The model- 
independent higgs'-strahlung process (see Figure [2jb)) 
where a scalar LLP is directly produced, and the signif- 
icantly larger resonant process shown in Figure |2|a), in 
which the A' can decay into LLP's. Data from past ex- 
periments such as the axion search at CHARM [16 have 
not been interpreted in terms of these processes anywhere 
in the literature. 

We note that electron beam-dump experiments like the 
SLAG experiments E137 [17] and E56 [H [19] are not 
competitive because of their lower energy and resulting 
decreased acceptance [53]. We also note that Hd can 
be directly produced from pion decay, or through pion 
decays into A' s followed by decays into ho- These mech- 
anisms are important in proton beam dumps where pion 
production is tremendously large for masses beneath 7X1^^. 
We refer the reader to [T0| for further analysis. 



A. Limits from CHARM Axion Search 

The CHARM axion search [16] had a 4 x 4 m^ detection 
area, extending over a 35 m decay region, 2.5 m deep 
and located 65 m behind a shielded proton target. A 
veto counter was placed in front of the detection volume, 
so that the search was sensitive only to exotic-particle 
decays within the detector volume. In all, = 2.4 x 
10^^ protons at an energy of 400 GeV were dumped on a 
copper target. No candidate e+e~, or 77 events 

were observed. 

We estimate event rates in the approximation that they 
are dominated by collisions in the first nuclear interaction 
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length X^^^ of the material: 

AT vnuc. 

Ax Ig/mol ^ 

= crn{X) X 2.0 X 10^pb"\ (2) 

where crn{X) is the per-nucleon cross-section for the pro- 
cess of interest, A = 63.5 is the average atomic mass of 
Copper and A^o is Avogadro's constant. 

We have modeled production in MadGraph [20] and 
used a Monte Carlo to determine geometric and lifetime 
acceptances. In the mass range 0.5 GeV < ttia' < 4 
GeV, the cross-section for resonant A' production for a 
400 GeV incident proton is well approximated by 



<T(m^O -20pb(^) (- 
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(for masses above 4 GeV, (j{mA') ~< cr(4GeV) x 
(m^//4 GeV)^). In addition, the angular acceptance 
varies considerably with the A^ mass. The typical en- 
ergies of A' produced resonantly were 40 — 50 GeV 
(with 10-20% exceeding 100 GeV); this is sufficient for 
a significant fraction of decay products to travel in the 
direction of the detector. Taking into account angular 
acceptance and decay probability. Fig. [3] shows the re- 
sulting limits on e^Br{X for i = e,/i as a 
function of LLP lifetime for various A^ masses. These 
limits were obtained assuming a decay A' hodD with 
^/iD — — -4771^1/ , but are not very sensitive to the 
detailed decay mode away from mass thresholds. 
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FIG. 3: Limits on e^Br{X as a function of decay 

length assuming resonant production of a vector A' , with sub- 
sequent decay to X. From bottom to top, curves correspond 
to ten events expected (none were observed) for uia' =0.6, 
1, 2, 3, 4 GeV and mhjj = OAttia'- The regions above the 
curves are excluded. The gray band represents the expected 
band for the scalar model of Section [lV| with the upper and 
lower lines corresponding to (ttia' jTUhjj) = (0.6,0.24) GeV 
and (4, 1.6) GeV, respectively. The width of the band comes 
predominantly from the opening of kaon decay channel at 
higher masses. 

For the higgs'-strahlung reaction, the cross-section is 
smaller by 4-5 orders of magnitude due to the combined 



effects of a weak coupling, phase space, and decreasing 
parton luminosity. For our benchmark rrih^ = 0.4mA' 
and 0.5 < m^/ < 10 GeV, the cross-section can be 
parametrized as 

aihoA') ^ (0.001 pb)(m/GeV)-^ (l + (m/4GeV)^)"^ 

(4) 

for = aEM' Limits on e^Br{X £~^£ ) for this 
reaction are shown in Figure [4] 

T(S) 




FIG. 4: Limits on e^Br{X £'^£~) as a function of decay 
length, assuming radiation of X off an off-shell as in the 
higgs^-strahlung process of 7 . From bottom to top, curves 
correspond to ten events expected (none were observed) for 
tua' = 0.6, 1, 2, 3, 4 GeV and ruh^ = 0.4mA/. The regions 
above the curves are excluded. The gray band represents the 
expected band for the scalar model of SectionjlA] with the up- 
per and lower lines corresponding to (m^/ , '^hjj) — (0.6, 0.24) 
GeV and (4, 1.6) GeV respectively. The width of the band 
comes predominantly from the opening of kaon decay channel 
at higher masses. 



III. DARK MATTER ANNIHILATION IN THE 
SUN AND EARTH 

We now move on to consider DM capture and annihi- 
lations in the Sun [22] and the Earth [23l [24], and 
the corresponding limits on the production rate and life- 
times of LLPs. We begin by reviewing the capture pro- 
cess of DM in celestial objects and the expected anni- 
hilation rates in relation to model parameters. Direct 
annihilation into SM particles can only be detected by 
looking for high energy neutrinos coming from the Sun 
and such searches have been discussed in the literature 
extensively. Annihilations into long-lived neutral parti- 
cles, however, may have very different signatures which 
we now consider. First, they may again be observed in 
neutrino detectors if the LLPs decay into muons for ex- 
ample. Second, they may be observed in electron de- 
tectors such as FERMI. Third, gamma-ray observations, 
such as EGRET, Milagro and FERMI can detect the ra- 
diation associated with electronic decays of the LLPs or 
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their direct decay into photons. This wide range of pos- 
sible probes of LLP decays is particularly welcome since 
the precise nature of these particles, if they exist, is still 
largely unconstrained. 



section below the experimental bound we find. 
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A. Production Rate from Capture and 
Annihilation of DM in the Sun and the Earth 

We consider several mechanisms for Solar and Earth 
capture and use the model described by Eq. ([T]) as a 
benchmark. We first describe the well studied case of 
elastic scattering of DM against matter and the associ- 
ated capture rate. We then examine inelastic scattering 
and the corresponding modified capture rate. In equilib- 
rium, the annihilation rate of DM is directly related to 
the capture rate. However, the fraction of annihilation 
into LLPs is model-dependent and so we parametrize it. 

The capture rate depends on the scattering cross- 
section of DM against nuclei. In the models we consider, 
the scattering is mediated through A' exchange and is 
given by [2,, 
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where Z and A are the atomic and mass number, respec- 
tively. Since the elastic scattering of DM against nuclei 
is constrained to be < 10~^^cm^ by direct detection ex- 
periments [25, 26 , the mixing is bounded to be e < 10~^ 
for rriA' ~ GeV. These bounds are dramatically weak- 
ened if DM predominantly scatters inelastically off nu- 
clei as in [27 . For example, DM-nucleon cross-sections 
cF^n ^ 10~^^cm^ can be compatible with direct detection 
limits and the DAMA modulation signal for splittings 
- 100 keV [28,. 



Elastic Capture 

Elastic capture in celestial objects proceeds through 
the basic mechanism discussed in [29, 30 . The scatter- 
ing can be spin-independent and spin-dependent. Both 
types of processes are constrained by direct-detection ex- 
periments, but the latter much less so. 

If spin-independent scattering cross-sections saturate 
bounds from direct-detection experiments [25l |26], the 
capture rates are 
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in the Sun and Earth, respectively [31]. The spin- 
dependent cross-section for proton coupling is most strin- 
gently constrained by KIMS [32 and setting the cross- 



In e /as tzc Capture 

Inelastic DM (iDM) was proposed in [27 as an expla- 
nation for the discrepancy between the annular modu- 
lation signal reported by the DAMA collaboration and 
the null results reported by other experiments. The ba- 
sic premise of this scenario is that dark matter scatters 
dominantly off nuclei into an excited state with splittings, 
~ 100 keV, similar to the kinetic energy of DM in the 
halo. With this hypothesis, the DAMA data is nicely fit 
with WIMP-niic/eon cross-sections of cr^x ^ 10~^^ cm^ 
(see [28 for details). While interactions with matter in 
the Earth are kinematically forbidden, the Sun's grav- 
ity provides the necessary kinetic energy to overcome the 
inelastic threshold and allow for capture of DM in the 
Sun. 

The capture rate in the Sun is shown in Fig. |5] for 
5 = 125 keV for several values of the WIMP-nucleon 
cross-section (see [33^ for more details). The parameter 
space is rather limited because above 5 ~ 500 keV there 
is no element in the Sun against which DM can scatter 
inelastically. We note that the high capture rates shown 
in Fig. [s] as compared with Eq. (|6| are mostly just a 
result of the higher WIMP-nucleon cross-section allowed 
in iDM models. 
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FIG. 5: The inelastic capture rate of DM in the Sun against 
the DM mass. The black (solid) curve correspond to an in- 
elastic model with 5 — 125 keV and (Tnx — 10"*^^ cm^, the 

10~^° cm^, and the purple 



blue (dashed) curve to cTnx - 
(dash- dot) curve to cfnx = 10" 
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Annihilation Rate 



detected near the Earth is, 



In equihbrium the annihilation rate of DM is exactly 
half its capture rate. More generally it is given by [34 , 



rA = -Ctanh' (t/TeJ, 



(9) 



where Teq is the equilibrium time and t ^ lO^^s is the dy- 
namical time of the system. In the elastic scattering case, 
DM can typically scatter in the Sun on a much shorter 
timescale than it annihilates. Hence, it thermalizes with 
the rest of the matter in the Sun and concentrates in the 
inner core as it approaches its equilibrium configuration. 
This is not the case for the Earth where t/r^^ <C 1 and the 
annihilation rate is suppressed with respect to the cap- 
ture rate by 2-3 orders of magnitude. If the annihilation 
cross- section enjoys an enhancement due to the Sommer- 
feld effect then this suppression can be tamed [35] . 

The situation is quite different for iDM since no more 
than a few inelastic scatterings can take place before it 
become energetically impossible to scatter any further. If 
the elastic cross-section is greater than cr^ ^ 10~^^ cm^, 
there are enough elastic collisions to ensure that equi- 
librium is reached. However, if <C 10~^^ cm^, the 
collision rate becomes so low that the WIMPs do not 
thermalize in the Sun. In Ref. [33 it was shown that 
equilibrium is nevertheless obtained, and we will there- 
fore equate the annihilation rate with half the capture 
rate in this case as well. 



B. Limits from the Sun and the Earth 

Bounds from Neutrino Detectors 

If the LLP decays into a pair anywhere between 

the Sun's surface and the Earth, these muons can be 
observed in underground neutrino detectors (Super-K, 
BAKSAL) as well as the neutrino telescopes (Ice-Cube, 
Ant ares) in either of two ways. If the decay happens 
very close to the detector (< km) the muon pair may 
be observed directly through their Cherenkov radiation 
as nearby tracks. Alternatively, if the LLP decays into 
a pair far from the detector then the muons will 

yield high-energy neutrinos which can be detected as they 
convert back into a muon in the rock or ice near the 
detector. 

The strongest bounds on the fiux of upward go- 
ing muons come from Super-Kamiokande and are 
O(IO^) km-2 yr-^ [36 . Fig. [g] depicts the bound on 
the annihilation rate of DM into LLPs that subsequently 
decay into muons, and the expected improvement in sen- 
sitivity with Ice-Cube and Antares. Some LLP's scatter 
or decay inside the Sun, while others do not decay before 
reaching Earth, or decay into modes that do not pro- 
duce neutrinos. Thus, the fraction of decays that can be 



fdecay = Br(LiP ^ /x's)e-^o/^' 



(10) 



where 7 is the LLP's boost factor, and 
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is the LLP mean free path in the Sun. The suppression 
^-Re/iint ~ I except for high e and low m^/ , to which the 
CHARM analysis is sensitive. When the lifetime becomes 
much smaller than the Sun's radius, the muons never 
escape and no useful limits can be derived. 

LLPs produced in the Earth must decay near (jl km) 
the detector, or else the muons are stopped in the Earth. 
In this case, the muons can be observed directly, although 
one must take into account the special nature of such 
events. Since the LLPs are very boosted, their decay 
products (namely the muons) are highly collimated with 
an opening angle of about ~ rriLLp/m^. The muons are 
therefore no more than about a meter apart when they 
reach the detector. We leave it for future work to inves- 
tigate the efficiency of observing such muon pairs, but 
there is no obvious reason to suspect that they should be 
considerably more difficult to see than ordinary events. 
Hence, the limits from Super-Kamiokande can be used to 
constrain the annihilation rate of DM in the Earth's core 
as shown in Fig. [?) Future results from Ice-Cube may im- 
prove on these limits by about two orders of magnitude. 
We note that such improvements will allow capture rates 
in the Earth (Eq. ^) to be probed despite the out-of- 
equilibrium suppression in Eq. (I9|). 
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FIG. 6: The purple curves show the current bounds from 
Super-K on the annihilation rate of DM in the Sun into LLPs 
assuming 100% branching ratio into 1-step (solid) or 2-step 
(dashed) cascades which result in muons. The red (orange) 
curves are the expected bounds from Antares (Ice-Cube). 
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FIG. 7: The purple curves show the constraints from Super- 
K on the annihilation rate of LLP in the Earth as a function 
of their lifetime for — TeV (solid) and = 0.3 TeV 
(dashed) . The orange curves show the expected improvement 
in sensitivity with Ice-Cube . The mass dependency reflects 
the improvement for higher masses, but does not take into ac- 
count possible degradation of the signal due to the colinearity 
of the muon pair. 



7 Ray Constraints 

Electronic decays of the LLP will usually lead to final 
state radiation (FSR) of photons unless the LLP is right 
above the electron pair threshold. This radiation can 
be detected by observations of the 7-ray spectrum of the 
Sun. Both EGRET and FERMI have looked at the Sun's 
emission spectrum for > 100 MeV and these agree 
well with the predicted spectrum from inverse Compton 
scattering of the Sun's light against cosmic rays 1401 HI] • 
We can therefore derive bounds on the electronic de- 
cays of LLPs by demanding that the resulting FSR is 
smaller than the observed gamma ray spectrum. As in 
the eletronic flux, the FSR energy spectrum is harder 
than that observed and so the highest energy measure- 
ments {Ej > 1 GeV) actually provide the most stringent 
limits. For that purpose, the measurement of the very 
high energy photon spectrum made by Milagro are also 
useful [42 . To obtain a bound on the flux we require 
the expected number of events in Milagro to be less than 
4791. The resulting bounds on the annihilation rate from 
the different observations are shown in Fig. [Sj 



Bounds from Electron/Positron Detectors 

Injection of highly energetic leptons between the Sun 
and Earth can also be observed in dedicated satellite 
missions that look for electrons/positrons in cosmic rays 
such as PAMELA [37] and FERMTLAT [38^. The en- 
ergy spectrum of such electrons is determined by their 
production mechanism only and does not suffer from the 
uncertainties present in galactic propagation models. We 
consider 1-step and 2-step processes, where the LLP's 
decay directly into electron/positron pairs or do so via 
an additional intermediate state, respectively. To obtain 
the energy spectrum for each type of cascade we assumed 
the different mass scales involved in the cascade are well 
separated and that all particles decay isotropically. The 
relevant formulae can be found in the appendix of [39] . 
Decay of the LLP into a pair of muons is very similar and 
is in fact well approximated as a 2-step process. 

For such short cascades, the electron energy spectrum 
is much harder than the cosmic-ray background. There- 
fore, the strongest constraints are derived from the high- 
est energy bins available. It is then straightforward to 
derive the bounds on the annihilation rate in the Sun 
by demanding that the resulting flux is lower than the 
highest energy bins in FERMI/LAT, 
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for the 1-step and 2-step cascades, respectively. We note 
that these are about an order of magnitude stronger than 
the neutrino constraints from Super-K. 



10 

^ 0.1 

X 

> 0.001 
o 

^ 10-5 
10-" 



100 200 



500 1000 2000 
M^(GeV) 



5000 ixlO^ 



FIG. 8: FSR constraints on Fann x fdecay from Milagro 
(green) , EGRET (purple) and Fermi solar data (blue) for one- 
step (solid) and two-step (dashed) cascade decays, with decay 
efficiency fdecay defined in (10). The bounds were derived as 
explained in the text. 



IV. EXAMPLE: DARK-SECTOR SCALARS 

In this Section we discuss the implications of the previ- 
ously identified limits for a particular class of long-lived 
particles: Higgs-like dark scalars that decay through 
loop diagrams with a rate controlled by the gauge ki- 
netic mixing e. We then discuss the dependence of these 
contours on variation of model assumptions. 

In a single-Higgs model the width of the light scalar 
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hf) to each lepton species is, 



rrih, (14) 



where I2 is a loop integral defined in [7] , and very nearly 
equal to 2/?^, where f3 = (l-4m|/m^^)^/^. This formula 
must be corrected by mixing angles in a model with two 
or more Higgses in the dark sector, but given our com- 
plete ignorance of the dark sector's structure we ignore 
these and take ([l4| as a benchmark. For decays above 
the muon threshold, this gives 

cr« (2x10^ cm)x ^^^,^^^,^_) ^ 

(^)- (15) 

We take the ratio of hadronic to muonic widths from a 
spectator model, as in [43j (we use the parameter choices 
of [4T with r = 1) - the branching fraction varies 

between 20 and 50% below the muon threshold, and falls 
to 1-2% above the kaon threshold. 

The terrestrial direct production limits and those from 
DM annihilation in the Sun (assuming annihilation prod- 
ucts include the A' and/or hn) are summarized in Figure 
[9j where for definiteness we have assumed rrih^ = rria^ = 
OArriA'^ and ajj = a. 

The most constraining terrestrial experiment is the 
CHARM beam dump, which searched for e+e~ and 
events and found no candidates in 2.4 x 10^^ pro- 
tons dumped. These results imply that e > 10~^ are ex- 
cluded for gauge boson masses up to a few GeV. Larger 
e ~ 0.1 — 1 is constrained by other measurements, specif- 
ically the muon g — 2. 

If weak-scale dark matter (m =100 GeV-1 TeV) an- 
nihilates through the A' ^ smaller e can be probed by 
measurements of the solar electron, photon, and neu- 
trino fiux. For TUhD < 2mr-, the photon fiux from final 
state radiation of dark Higgses that escape the Sun is 
most constraining, with sensitivity to annihilation rates 
~ 10^^ — 10^^ /s, four to five orders of magnitude below 



the expected capture rate for a 
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assume the capture rate is determined by the scattering 
cross- section (|5| from exchange of an A' ^ with no mixing- 
angle suppression. In this case, the solar production rate, 
like the terrestrial beam-dump production rate, scales as 
e^, leading to a very similar shape of exclusion region 
(blue region for tudm = 100 GeV, dark blue dashed line 
for rriDM = 1 TeV). 

Unlike muons, pions, and kaons, the r decays 
promptly, producing neutrinos that can escape the Sun 
even if the dark-sector Higgs decays promptly. Thus 
dark matter scenarios above this mass threshold are con- 
strained (green vertical band in Fig. |9|, with little sen- 
sitivity to dark matter mass. 

As we have noted, the decay lifetimes of Higgs-like 
scalars are quite model-dependent, and can be signifi- 
cantly longer or shorter than those quoted. This would 
qualitatively affect both beam-dump (red) and solar FSR 
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FIG. 9: Combined limits on a hig gs-li ke scalar hn produced 
in A' decays, with leptonic width (14), m^/ = 2.5m^^, and 
ttD = tt. The gray region in the upper left is the model- 
independent limit on e for given A' mass from (g — 2)^ (see 
[45]). The red region corresponds to the limit from direct 
A' production in CHARM, with ho decays into electrons or 
muons (see Sec|ll|. Where this limit is obscured by others, 
the outline of the CHARM exclusion is drawn in dark red. 
The blue and green regions are solar limits (Sec. |III B|) rele- 



vant if weak-scale dark matter couples to the A\ and kinetic 
mixing mediates its scattering off ordinary matter. The blue 
and green regions are excluded by photon and neutrino flux 
observations, respectively, for a benchmark dark matter mass 
= 100 GeV. The dashed blue line outlines the analogous 
exclusion from photon flux when rrix — 1 TeV. The neutrino 
limits are controlled by the r threshold and are insensitive 
to the dark matter mass. Possible strategies for exploring 
regions I, II, and III are discussed in Section [y] 



(blue) constraints. The upper limits on e correspond to 
typical 7cr comparable to the thickness of the beam- 
dump shielding (or Rq)^ with only logarithmic sensi- 
tivity to production rate. These shift up (down) for 
longer (shorter) lifetimes. The lower limits are deter- 
mined equally by the lifetime (which typically exceeds 
the size of the experiment) and production rate, so it is 
somewhat less affected by variation of the decay length. 

The solar constraints also depend on various assump- 
tions about dark matter. The dependence on dark mat- 
ter mass, shown in the Fig. [9] (dashed blue line), arises 
mostly from the resulting boost of the late-decaying Hig- 
gses. In models with a reduced DM-matter scattering 
cross-section, the luminosity of the Sun in dark matter 
annihilation is reduced and the resulting FSR limits are 
weakened. We note that a mixing angle weakens the 
lower-e limit significantly, but has only logarithmic im- 
pact on the high-e sensitivity limit. 
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V. DISCUSSION 

Dark sectors coupled to the Standard Model through 
gauge kinetic mixing may contain an array of macroscop- 
ically long-lived particles (LLP's) in addition to the ki- 
netically mixed gauge boson. Gauge bosons can escape 
detection in accelerator experiments if they decay into 
LLP's with lifetimes between 10 cm and 10^^ cm, rather 
than to Standard Model fermions. LLPs with lifetimes 
in this range would escape collider detectors unobserved, 
but are sufficiently short-lived that BBN constraints are 
largely irrelevant. 

Remarkably, two classes of data constrain these sce- 
narios across these fourteen decades in lifetime. The pro- 
ton beam dump experiment CHARM places strong con- 
straints on scalar LLP's with 10 < cr < 10^ cm decaying 
into muons or electrons. If dark matter annihilates into 
dark sector LLP's, then Solar capture of dark matter re- 
sults in limits on LLP's with 10^ cm ^ cr < 10^^ cm for 
similar decay modes. 

These limits each depend on different assumptions and 
production mechanisms, so they cannot be combined into 
a robust exclusion. Nonetheless, the case of a Higgs- 
like boson in the dark sector decaying through kinetic 
mixing allows an illustrative comparison. The sensitiv- 
ity of solar limits — when applicable — nearly matches 
on to the parameter region excluded by CHARM. To- 
gether, these constraints disfavor A' decays to scalar s, 
pseudoscalars and other long-lived particles in the param- 
eter region best suited to explain the DAMA modulation 
signal through iDM [Ml [151 SS]- In particular, if LLP's 
decay through higher-dimension operators, lowering the 
mass scale of these operators to evade BBN constraints 
may bring them into conflict with experimental data. 

Given these limits and the remaining unconstrained 
parameter regions, our analysis suggests several promis- 
ing searches with existing data and small-scale experi- 
ments: 

Prompt Vector Decays to £~^£~: Remarkably, for 
e ~ 10~^ - 10~^, A' decays into LLPs are 
much more constrained than decays directly 
into Standard Model states. This observation 
compounds the motivation for a sharply defined 
search program to constrain direct A' decays to 

£^£- pig^. 

Searches in e+e~ Collider Data: LLP lifetimes too 
short to be observed in a beam-dump experiment 
{jCT < 10m at Ellp ~ 50GeV, e.g. region I 
of Fig. [9| give rise to observable displaced decays 
within the detectors of e+e~ collider experiments 
such as Babar, Belle, KLOE, and BES-III. More- 
over, such lifetimes arise from e > 10 ~^ for which 
many thousands of events are possible in j -\- A' 



and hnA' production. The most visible such decay 
may be into (or tt+tt") pairs reconstructing 

a displaced vertex, which may be accompanied by 
a photon, invisible particles, and/or additional e, 
/i, or TT pairs. 

High-Energy and High Intensity Beam Dumps: 

High-energy beam dumps with cumulative charges 
on target exceeding a coulomb may extend sensitiv- 
ity into region II of Fig. [9j Neutrino factories with 
near detectors such as MINOS have geometries 
well-suited for beam-dump LLP searches |47] , 
and integrated luminosity ~ 10^ larger than that 
of CHARM. 

Solar Gamma Ray and Electron Measurements: 

LLP lifetimes ~ 10^ — 10^ seconds (region HI in 
Fig. |9| are too short to disrupt light element 
abundances, but only a fraction of LLP's can 
decay between the Sun and the Earth, so very 
sensitive solar measurements are needed to probe 
these scenarios. Solar gamma-ray observations 
above 10 GeV from ACTs and FERMI can explore 
this region, as can searches for electronic signals 
correlated with the Sun [48] . 

Milli- Charged Particle Searches: Electron beam- 
dump experiments searching for milli-charged 
particles (e.g. [49l |50]) may be sensitive to LLP 
scattering off a nuclear target in the active 
detector volume, such as ho -\- N A^ -\- N. 
Searches in existing data may be sensitive to 
e > 10-^(100 MeV/mA') and tullp < 100 MeV. 

Note added: While this work was brought to completion 
two related works [51], and [52] appeared. 
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